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We study the statics and the dynamics of domain patterns in proper hexagonal-orthorhombic fer- 
roelastics; these patterns are of particular interest because they provide a rare physical realization 
of disclinations in crystals. Both our static and dynamical theories are based entirely on classi- 
cal, nonlinear elasticity theory; we use the minimal theory consistent with stability, symmetry and 
ability to explain qualitatively the observed patterns. After scaling, the only parameters of the 
static theory are a temperature variable and a stiffness variable. For moderate to large stiffness, our 
static results show nested stars, unnested stars, fans and other nodes, triangular and trapezoidal 
regions of trapped hexagonal phase, etc observed in electron microscopy of Ta4N and Mg-Cd alloys, 
and also in lead orthovanadate (which is trigonal- monoclinic) ; we even find imperfections in some 
nodes, like those observed. For small stiffness, we find patterns like those observed in the mineral 
Mg-cordierite. Our dynamical studies of growth and relaxation show the formation of these static 
patterns, and also transitory structures such as 12-armed bursts, streamers and striations which 
are also seen experimentally. The major aspects of the growth-relaxation process are quite unlike 
those in systems with conventional order parameters, for it is inherently nonlocal; for example, the 
changes from one snapshot to the next are not predictable by inspection. 

PACS numbers: 81.30.Kf, 68.35.-p, 62.20.Dc, 61.70.Ng 



I. INTRODUCTION 



Ferroelastica^l'cl are solids that undergo diffusionlcss, 
shape-changing phase transformations, usually first- 
order. The low-temperature (low-T) phase is, in prin- 
ciple, one of several degenerate products ("variants") 
differing only in orientation. But one rarely sees the 
homogeneous product phase; instead, multiple nucle- 
ation events, constraints, intrinsic inhomogeneities, etc, 
produce domain structures composed of multiple vari- 
ants separated by domain walls (ideally with certain 
orientations). 

This article deals with hexagonal-orthorhombic (H- 
O) ferroelastics, where the three equivalent stretch di- 
rections give rise to three variants. Examples are the 
mineral Mg-cordierite (Mg2Al4Si50i8, References 
Ta4N (Reference 0), and Mg-Cd alloys (References] 
The most interesting feature of these materials is that 
they pravide a rare physical realization of disclinations in 
crystala3o. The disclinations result because domain walls 
rotate the variants as well as link them, all nodes but one 
failing to close to 360° ; the consequence is the fascinating 
domain patterns observed in these materials and in the 
closely related trigonal-monoclinic material lead ortho- 
vanadate [Pb3(V04)2, References p^ , pd| . These patterns 
differ strikingly from those in cubic-tetragonal (C-T) and 
tetragonal-orthorhombic (T-0) systems; for example, the 
star disclination, often self-similarly nested, has no coun- 
terpart in the other systeHi&r-, 

Of the several theoriegl3"t^ dealing with these materi- 
als, only that of Wen et alr3 found patterns consistent 
with some of those observed. Their phase-field kinetic 
theory is however highly phenomenological, being based 



on order parameters with no direct physical interpreta- 
tion; an expansion to sixth order is necessary, and even 
the scaled theory requires at least five parameters. 

We investigate H-0 ferroelastics by expanding the 
energy density in physical quantities only, namely the 
strains and their derivatives, to fourth order in the pri- 
mary order parameter. Computational resources allow us 
to consider only two-dimensional structures, but this is 
not much of a restriction, for these are what the literature 
reports. Our density respects the symmetry of the high-T 
phase; it is minimal, for it contains only essential terms. 
Scaling reduces the parameters in the density to two, a 
temperature variable and a stiffness variable, both read- 
ily determined from experiment. Our calculated struc- 
tures depend strongly on the stiffness, which determines 
how the system responds to the disclinations generated 
by wall junctions. We investigate both the statics and 
the dynamics of the theory, using in both methods which 
satisfy the compatibility relations exactly. 

Our static theory, described in Section II, is based 
solely on classical elasticity theory; it includes nonlin- 
ear and gradient terms but no additional phenomenol- 
ogy. Section III solves for the two kinds of wall linking a 
given pair of variants. Section IV describes static struc- 
tures found numerically. Section V first describes the 
equations of motion of our dynamical theory (also based 
solely on classical elasticity theory) and then goes on to 
present structures as they evolve in time. 

Although our theory is strictly valid only for proper 
ferroelastics, where the strain is the primary order pa- 
rameter, we find that it is sufficient to explain qualita- 
tively the patterns observed in all four materials listed 
above, namely several kinds of three-fold nodes, a com- 



1 



mon four-fold node, stars (both nested and unnested), 
the fan node and another untitled node, and even im- 
perfections in some higher-order nodes. Near Tc, we find 
triangular and t rape zoidal pockets of the high-T phase, 
also as observedll^llli. Many of these structures were not 
found by Wen et al£3. Our dynamical studies start from 
the perturbed high-T phase. In the growth stage, we see 
12-armed bursts as observed in Mg-cordierite. In the re- 
laxation stage, but still far from the converged state, we 
see streamers and striations also as observed, and also 
the formation of the structures seen in our static stud- 
ies. These are the first simulations to report the bursts, 
streamers and striations. 



II. STRAINS AND STRAIN ENERGY 

Hexagonal-orthorhombic ferroelastics require five 
strains in general. Since we consider only columnar struc- 
tures uniform in the 3 direction (normal to the basal 
plane), we need only the dilatational strain ei and the 
deviatoric/shear strains 62 and 63. In the linear approx- 
imation for the strain tensor. 



ei = (wi,i + ^2,2) /2 

62 = (-"1,1 - "2,2) /2 

63 = ('"1,2 + "2,1) /2 



(la) 
(lb) 
(Ic) 



where u(cci, X2) is the displacement of the material point 
originally at x, and Uij — djUi = dut/dxj. 

The free energy F, per unit length in the 3 direction, 
is the integral F = J^T cPx oi the free-energy density 
J- over the undeformed area A. Since the following is 
the first study of structures in proper H-0 ferroelastics, 
and since there is no evidence that other terms are re- 
quired, we use the simplest possible density consistent 
with symmetry, stability and ability to explain qualita- 
tively the structures observed. Of course the terms omit- 
ted may be numerically important; in fact a class of these 
terms is qualitatively important, for we show below that 
it breaks an important degeneracy. 

The minimal— density is slightly simplified from one 
used previouslyllj: 



B2 
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(2) 



All coefficients but A2 are positive and independent of T. 
The Ai and A2 terms are the contributions from homoge- 
neous volume changes and shears respectively, the latter 
in the high-T phase only; the parameters Ai and A2 are 
related to the Voigt constants by Ai = 2(Cii + C12) 
and ^2 = 2(Cii - C12) = iCee. But A2 softens with 
decreasing T, as A2 — a{T — Tg), where Tg is the stabil- 
ity limit of the H phase and a is a material-dependent 
constant. The instability due to the vanishing of A2 is 



isotropic. The orientations of the three O variants are 
selected by the anisotropic B2 term, as described below. 
The quartic C2 term, isotropic, provides stability. Since 
ei is not the primary order parameter, the term ^Aie\ 
is often omittedtHHlil. We show below however that this 
term plays a crucial role in determining the morphology; 
very different structures are found for large vs. small Ai . 
The d2 term, which penalizes spatial inhomogeneities in 
the primary order parameter, is necessary to give the 
domain walls an energy and so prevent subdivision into 
arbitrarily fine variants. A discussion of other possible 
derivative invariants is given in the Appendix. 

The H-0 transition occurs at A2 = |i3|/C2; it is first- 
order. For A2 > JB2/C2, the energy has only the H min- 
imum, with 61 = 62 = 63 = 0; this minimum disappears 
for A2 < 0. For A2 < JB2/C2, there are in addition 
three degenerate O minima with 61 = 0, symmetrically 
located in the (62,63) plane: 



where 620 = 



62 = 620 , 63 = 

62 = -620/2 , 63 = -\/362o/2 
62 = -620/2 , 63 = -t-\/362o/2 

1/2 



(3a) 
(3b) 
(3c) 



B2 



BI-AA2C2) ' /(2C2) ; at T„ 



620 — I-B2/C2- The three variants are elongated in the 
i, — ii -I- ^2, and — ii — ^2 directions respectively. 



III. VARIANT- VARIANT WALLS 

We seek the solution linking the two variants of Eqs. 
(3b) and (3c), the results for this pair being simpler than 
for the others; the analysis is easily extended to treat a 
band of walls. The key new results are the expressions 
(7) for 61 and W3 in terms of 62 and 63. 

We start from the first-order compatibility relations 
Wi,i2 — Mi, 21; with the linearized strain tensor, these are 



82 (61 -t- 62) = di (63 +UJS.) 
di (ei - 62) = 82 (63 - 0^3) 



(4) 



where W3 = (ui^2 — 'U2.i)/2 is the z component of the local 
rotation a;. Elimination of gives the more familiar 
second-order relation 

{df + dl) 61 - {dl - dl) 62 = 291^263 . (5) 

On applying the boundary conditions 

e2(±oo) = -620/2 , 63(±cx)) = ±V362o/2 , (6a) 



6i(±oo) = , Ci;3(±oo) — ±rt , 



(6b) 



and demanding that the strains and L03 be functions 
of the single variable X = xi cos/3 + X2 sin/3, one finds 
easily the two solutions 

/3 = , X = xi , 61 = +(62 + 620/2) , CJ3 = -63 (7a) 
/3 = 7r/2 , X = a;2 , 61 = -(62 + 620/2) , UJ3 = +63 (7b) 
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Then, with the other two pairings, one has a total of 
six possible wall orientations, at multiples of 7r/6 to the 
hexagonal 100 axig3. These results are independent of the 
details of the free-energy density; we made only three as- 
sumptions: the strains are small, they (and lo^) are func- 
tions of a single variable, and they have the asymptotic 
values (6). 

The|-|two kinds of wall are orthogonal, as already 
knowiH. They are obviously different in nature; it is 
enough that the angles between the stretch axes are dif- 
ferent, but we note also that the dilatational strains have 
opposite sign in the wall region. The wallsjieed not pos- 
sess the twin property; in Mg-Cd alloysHtl, the "reflec- 
tion" walls are twin walls but the "rotation" walls are 
not. 

The rotations in Equations (7) are the key to under- 
standing the domain patterns. 

For (3 = 0, X = xi: the variant at Xi —>■ +oo, with 
63 > 0, is rotated negatively (ws < 0, clockwise); the 
other variant is rotated positively; the angle between the 
stretch directions is not n/3 but rather 7r/3 — ^/3e2Q- 
For P = 7r/2, X — X2'. the variant at X2 +oo, with 
63 > 0, is rotated positively (counterclockwise); ... ; 
the angle between the stretch directions is not 27r/3 but 
rather 2tt/3 — •\/3e2o- 

These rotations are well known in the metallurgical liter- 
ature, where they are obtained from geometrical consid- 
erations. In Mg-Cd alloys, the angles between the stretch 
axes of the linked jiariants are calculated from the lat- 
tice constants to ben 57.5° or 117.4°, in good agreement 
(± ^ 1°) with the observations^. The leading correc- 
tions to Eqs. (7) (due to nonlinear terms in the strain 
tensor) are easily found: 61 = ±(62 -f ^620) + O (6'^) and 
^3 = Te3(l T |e2o) +0 (e^). 

Our resultsJpr H-0 wails jare in some respects the same 
as for C-TBeI and T-oN^Ea systems. In all three sys- 
tems, a given pair of variants can be Jinked by one of two 
possible walls, at 7r/2 to each otherEl. A given wall not 
only links the two variants, but also rotates them differen- 
tially by an amount proportional to the order parameter, 
as seen for example in Eqs_-£7). Collisions of wall bands 
generate new kinds of wallscj, including a strange "wall" 
between two variants with identical strains but opposite 
rotations. In both C-T and T-0 systems, where the an- 
gle between the variant axes is tt/2 minus a correction, 
the two possible wall structures are identical with respect 
to the dilatational and shear strains; in fact, they are re- 
lated by a symmetry operation, and so have the same 
energy. In H-O systems on the other hand, there is no 
symmetry operation relating the two possible structures 
and so they generally have different energies. 

To determine the structure of the wall linking the vari- 
ants of Eqs. (3b) and (3c), we first use Eqs. (7) to 
eliminate 61 from Eq. (2) for the density: 





Here D2 — -D3 = ^2 for the density of Eq. (2), but D2 
becomes D2 — di + d2 if the di term is included in the 
density (see the Appendix); then the length scales of 62 
and 63 for the wall can in principle be different, even 
though their gradients appear with the same coefficient. 
Obviously the two kinds of wall {X = xi and X = X2) 
have the same energy. This degeneracy is however due 
to the particular form of Eq. (2); it is broken if ws-add, 
to T any invariant odd in 61, for example the termEirE2l 
Eei{e2 + 63). In principle, such terms can give one wall 
type a much larger energy than the other; if so, then the 
domain patterns may be quite unlike those in cordierite, 
Ta4N, Mg-Cd and lead orthovanadate. 

The compatibility relations are satisfied in Eq. (8), 
and so the Euler-Lagrange equations can be obtained by 
varying 62 and 63: 

Ai(62 -I- 620/2) + A2e2 - B2{el - el) + (7262(6^ -f el) 

= D2d^e2/dX'^ , (9a) 
A263 + 2B26263 + C263(62 + el) = Dsd'^es/dX^ , (9b) 

with the boundary conditions of Eq. (6a); the Ai term 
in Eq. (9a) results from the compatibility relations. An- 
other route to Eqs. (9) starts from requiring F, with the 
density of Eq. (2), to be stationary with respect to the 
displacement u; then the strains are assumed to depend 
only on xi or X2 and Eqs. (7) are used to eliminate 61; 
finally, the results are integrated to give Eqs. (9). 

Equations identical in appeasance to Eqs. (9) were de- 
rived and solved in our studylla of C-T ferroelas1ji«i. At 
A2 — — 2i3|/C2, they have an analytical solution£3 with 
62 = —620/2 and so 61 = 0; at any other T, however, the 
dilatational strain is not zero in the wall region. The an- 
alytical solution at A2 — — 2i3|/C2, numerical solutions 
at other T, and numerical and analytical results for the 
wall energy are given in Ref. |l^. These C-T results carry 
over to H-O systems with only changes in notation; we 
point out however that they seem to differ significantly 
from the results of Ref. |l^. 

IV. STATIC STRUCTURES 

To reduce the number of unknown parameters, we 
transform variables by 

ej [2 X 10^B2/{3C2)] e, , (10a) 

X, -> ^M2C2/{2Bl)x, , (10b) 

T ^[^y.lQ^Bl/{?>lCl)]T . (10c) 

Because we work with the displacement, the first 
(strain) scaling assumes small strains (so that the 
quadratic terms in the strain tensor can be neglected). 
We also define a dimensionless temperature by r = 



3 



(T — To)/{Tc — To), and a dimensionless anisotropy pa- 
rameter by C = Ai/A2{T = Tc) = 9AiC2/(2B|). In 
terms of the new variables, the density is 



c , . 



69 + 63 
1 

2 



- 3 (^2 - 3626^) 



(Ve2) V (Ve3 



(11) 



where 6 = 3 x 10^ and c ~ 2x 10^; the scaled theory then 
has only two parameters, r and C. The transition (T = 
Tc) is at r = 1, the order parameter at Tc is eaa = 10"'^ 
(an arbitrary value) and the analytical solutiontll is found 
at T = —9. To relate r to the physical temperature T, we 
require both Tc and Tq; the latter can be estimated from 
the T dependence of the order parameter (for example, 
the value at Tq is 3/2 the value at Tc). 

Wishing to study wall configurations, we forced walls 
into the systems by applying periodic boundary con- 
ditions to the displacement u; this apparently simple 
requirement is sufficient to generate many metastable 
states. Free boundary conditionai3, it seems, cannot 
generate the complex structures seen in experiment, and 
fixed boundary copditions (specifically u = on and out- 
side the boundaryn3'E2l) give results which suffer severely 
from finite-size effects. The static structures described 
in this section were found by minimizing the free energy 
with respect to the components of u, using a conjugate- 
gradient method; the iteration was continued until nu- 
merical errors prevented further decrease in the energy. A 
variety of starting configurations was used, often random 
displacements. Finite-size effects are certainly present 
in the systems we can study (they favour simpler struc- 
tures); we see the basic structures observed experimen- 
tally, however, and so these effects seem unimportant 
qualitatively. 

Detailed comparison with experiment is unfortunately 
not possible for any of the materials; in fact, experiment 
offers no guidance at all regarding the parameters r and C, 
of even the minimal theory. Accordingly, we investigated 
several sets of parameter values; the desired comparison 
can be made when the parameters are measured. 

As the temperature decreases, the order parameter 
(e2o) of the homogeneous (O) product phase increases 
and the wall width (proportional to l/e2o) decreases. 
The two curvatures at the O minima also change, but 
in different ways. These are 7|| and 7j_ in the directions 
parallel and transverse to the line from the origin to the 
O minima: 



7||(r) = r 

7_L(r) = T 



2^620 

26e2o 



-20 J 



(12a) 
(12b) 



both increase monotonically with decreasing T . For 
Tbk < T <Tc (-9 < r < 1), 7|| < jj_ with (7|[,7j.) = 
(1,9) at T = Tc. At T = Tbk (t = -9), the curvature 
is isotropic, with 7|| = 7_l = 27. For T < Tbk, 7|[ ^^i-^-L- 
The changing curvatures cause the wall trajectoriesta to 
start closer to the easy direction in the (62,63) plane. 



The point of introducing these quantities is that our 
structures depend strongly on the ratio of the curvature 
C (for dilatational strains) to the curvatures 7|| , 7j^ (for 
departures of the deviatoric/shear strains from values in 
a homogeneous O variant). This ratio is small for "soft" 
systems and large for "stiff" ones. This strong depen- 
dence on the relative stiffness arises, in the final analysis, 
because the walls rotate the variants as well as link them. 
Because of the rotations, disclinations are formed when 
walls collide to form junctions (or nodes). Since a sin- 
gle disclination cannot exist in a crystal, the walls form 
arrays of disclinations of both polarities, often in pairs 
(dipoles) and higher multipoles. Finally, soft systems 
and stiff systems respond very differently to the disclina- 
tions. In soft systems, the walls can deviate considerably 
from the optimal orientations found above; additional di- 
latational strains are then generated, but these have low 
energy since Ai is relatively small. In stiff systems, on 
the other hand, the walls are straight on the whole, and 
their directions are highly constrained; other means must 
be found to minimize the energy. 

To describe the nodes, we use the notation (fc, ?,•••) 
where the components fc, Z, • • • are successive angles in 
units of 7r/6 (only approximately); this notation, though 
incomplete because it does not distinguish between the 
two types of walls, is adequate for our purposes. The 
43 possible disclinations, including the star, have been 
catalogued and tbp values of the disclination angle.(5 given 
for Mg-Cd alloygJ; those with ^smallest |<5| foUowO: 

1. (5 = for the (1,3,5,3) nodetlil. 

2. \S\ — 0.243° for the star disclination (a three-pointed 
star with tip angle ^ 7r/6, often self-similarly nested), 
the (1,1,1,1,1,7) fan node, and the untitled (1,1,1,4,1,4) 
node. 

3. \6\ — 0.486° for the "starburst", our term for a node 
with 12 walls separated by ^ 7r/6. 

4. \S\ = 2.472° or 3.201° for the (1,1,1,1,1,1,1,1,4) node. 

5. \S\ = 2.715° or 2.958° for the (2,5,5), (1,1,10) and 
(1,4,7) three- fold nodes and for two five-fold nodes. 
Many of the 43 disclinations are not seen because \6\ 
is too large; for exaxople, the at-first-thought obvious 
(4,4,4) node fails bylS'tl ~ 8°. And one expects nodes 
with many arms to appear less frequently. 

Figure 1(a) shows a portion of a structure for a soft 
system {( = 10) at r = —50, where (7||,7_l) « (117,52). 
In this and other structures we found for the same pa- 
rameter values, we see only three-fold nodes and bends, 
even though (1,3,5,3), star, fan and other nodes have 
smaller disclination angles. In such soft systems, the 
walls can deviate markedly from the high-symmetry di- 
rections; they bend gently on the whole, but sharp bends 
can occur, especially at necks like that at the bottom of 
the figure. Examination of the rotation oj^ shows that 
the necking results from the collision of identicaJ^vari- 
ants with different orientations. As in T-0 systemsc2l, the 
necking reduces the length of the collision region where 
large dilatational strains are generated. At C = 1000, the 
neck and the 7r/2 bends were replaced by a (1,3,5,3) node. 
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though with other significant changes. Figure 1(a) shows 
also two well formed (2,5,5) nodes (prominent in Ref. ||); 
not shown are poorly formed (1,1,10) and (1,4,7) nodes 
(also seen in Rcf.J^|-.The overall morphology is like that 
seen in cordieritccTElB, but lack of detail in the published 
micrographs prevents a direct comparison; there seems 
to be room here for further experiment. 

Figures 1 (b) and 1 (c) show typical structures obtained 
also at r = —50, but with the moderately soft value 
C = 100. Nested and unnested stars appear, and also 
(1,3,5,3), (2,5,5), (1,1,10) and (1,4,7) nodes, though all 
are imperfect. Nested stars like that in part (b) are com- 
monly seen in experiment; the unnested star of part (c) 
is like that at the right of Figure 3(c) of Ref. |l^. 
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FIG. 1. Greyscale plots of {e'i + e'i)^^^ for selected re- 
gions of static structures. In the white regions, the devia- 
toric strains are near one of the three orthorhombic minima; 
the black regions are domain walls. Positions 2, 3, and 4 
show (2,5,5), (1,4,7) and (1,1,10) nodes respectively. Part (a) 
shows a soft system {(, — 10) at low temperature (r — —50). 
The neck at position 1 separates two domains with the same 
strains but different orientations. Parts (b) and (c), both 
moderately soft systems (^ = 100) at r = —50, show nested 
and unnested stars respectively. Part (d) shows a moderately 
soft system {( = 100) at a higher temperature (r = —20); 
the length scale is different. (1,1,1,4,1,4) and (1,3,5,3) nodes 
appear at positions 5 and 6. 

Figure 1(d) shows a (1,1,1,4,1,4) node, though the ver- 
tex in imperfect. This node, similarly imperfect (walls 
jog at the node), appears in Fig. 7(c) of Ref. ^ and in 
Figs. 3(b) and 3(c) of Ref. ^ We find that the node 
tends to split into simpler nodes at lower T, as the effec- 
tive stiffness (/j decreases. We often observed, but do 
not display, staircase structures and the "spiral branch 



configuration" shown in Fig. 3(a) of Ref. [To| . 

The structure of Figure 2(a), for a stiff system {( — 
1000) at T = —50, shows mostly objects with small discli- 
nation angles, namely (1,3,5,3) nodes and stars; the few 
three-fold nodes are badly distorted (1,4,7) nodes. 

Figure 2(b) shows a structure formed close to Tc. The 
most important features are the triangular and trape- 
zoidal pockets of the high-T (hexagonal) phase; the first 
is seen in Refs. |l^,|ll|, and the second in Fig. 1(e) of ReL 
pT| . Trapped high-T phase was not seen by Wen et al.tj 
who looked at a single temperature, perhaps too low. The 
pockets can be viewed as small grains trapped by rigid 
walls which prevent them from changing shape; they can 
transform to the low-T phase only through formation of 
internal walls, which however cost too much energy at 
this T. The structure of Figure 2(b) was generated by 
warming. In the initial structure, the stars were nested 
and the trapezoids also showed internal structure; both 
internal structures gradually faded away with warming. 
Figure 2 of Ref. |ri| shows comparable before-and-after 
pictures found on cooling. 

Figure 2(b) shows also the fan node, though the vertex 
is imperfect. Imperfect fan nodes are seen in Fig. 7(b) of 
Ref. ^ in Fig. 11(b) of Ref. || (the high-resolution image 
shows 3 three- fold junctions), and in Fig. 3(d) of Ref. |l^. 

We found structures also for ( = 0.5, 5 and 50, at 
T = 0, where (7||,7_l) = (4.5, 13.5); the patterns differed 
significantly in some respects from those described above. 
These results and those for r = — 50 suggest that the stiff- 
ness ratio C/7 is of order 1 for Mg-Cd at the temperature 
of the observations. 

Since C/7 depends on T, significant structural changes 
may result from changing T: if C/7 is large at Tc and 
small at low T, then stars and other structures charac- 
teristic of stiff materials may decay as T decreases. 



V. DYNAMICS 

We derive the equations of motion using the La- 
grangian formulation, starting from C — Sp^i J' where 
T — p(ui)^ 12 is the kinetic-energy densitytSEj. The re- 
sult is 



(13) 



where Gij — dC/dui^j is the stress tensor and a',^ is the 
dissipative stress tensor. The most general form for the 
latter, which describes non-conservative forces (friction) 
in the system, is u'^j = A^jf.iUk.i to lowest order; we note 
that the dissipation cannot depend on the velocity field 
ilk itself, only on its space derivatives. The dissipative 
stress tensor is conveniently found from the dissipative 
function by = d'if/duij. The form which respects 
the symmetry of the parent (hexagonal) phase is 



A' 



(14) 
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We assume that the dissipation is large and drop the 
inertial term pui] the equations of motion are then 

A[ei,i + A'^{e2s + 63.2) = -{Gi.i + G2.1 + €3:2) , 
A[ei,2 - ^2(62.2 - 634) = -(Gi,2 - G2,2 + G3,i) , (15) 

where Gi — SF/Set: 

Gi=^iei, (16) 
G2 = ^262 - B2{el - el) + C2e2{el + el) - d2'^^e2, (17) 

G3 = ^263 + 2526263 + G263(62 + 6^) - d2^'^e3. (18) 

ActuaUy, we have to use the equations of motion in terms 
of the displacement components: 

/ j^^dl + A'^ {dl + dl) A\d^d2 \( 
\ A',d^d2 + A'M + dl) J[u2j 

— —2 ( ^^'^ ^ ^^'^ ^ ^^'^ ) (19) 

\ Gi^2 ~ G2,2 + G34 J 

The time is scaled so that A'2 = 1; the dependence on 
the parameter A'-^ was found to be weak during the H-0 
transition, but the final states (completely relaxed con- 
figurations) did depend on A[. We used A'^ = for the 
results shown in Figure 3. 

These equations must be solved numerically. Each 
time step started with the displacement u(a;i, X2) known 
at the points of a 256^ grid. The space derivatives of Ui 
were approximated by finite-difference expressions cen- 
tered on a 5 X 5 grid. The non-linear terms on the right- 
hand sides were computed, and then the equations were 
fast-Fourier transformed to give algebraic equations. The 
solutions were advanced in time using the Euler method, 
and then fast-Fourier transformed back to position space, 
to start the next step. 

All our dynamical studies were started by perturb- 
ing the high-temperature state at a single point. The 
initial state certainly does not represent the nucleus of 
the transformation, but we expect its growth to resemble 
that in experiment; in partial confirmation, we see struc- 
tures much like those observed. The perturbation grows 
quickly with time, but eventually it collides with self- 
images generated by the boundary conditions; of course 
the collision stage also has no relation to experiment. 
The hexagonal (high-T) phase then disappears (unless 
the temperature is close to Tc); thereafter, we expect our 
simulations once again to resemble experiment (again we 
see structures much like those observed) . The relaxation 
stage which follows is quite slow. We observe the for- 
mation of the objects already seen in the static studies 
reported above, but more interestingly also well-defined 
structures (streamers, striations) seen infrequently or not 
at all statically, most of which eventually decay. 

In stiff systems (large C/i)j large domains can be pen- 
etrated by streamers (thin, wedge-shaped intrusions with 
relatively thick walls off the optimal directions); one of- 
ten sees also striated patterns of two alternating variants. 



The striations and streamers are relics of the H-0 trans- 
formation. Dilatational stresses occurring at the inter- 
face between the high-T and low-T phase arejplieved by 
local fine-scale twinning of the low-T phaseS. The re- 
sult is more complicated domain patterns (striations and 
streamers) shortly after the transformation is complete. 
Star embryos appear in the vicinity of triple junctions; 
these form near mis-oriented or curved domain walls after 
the surrounding area has transformed. In contrast, sys- 
tems with small C/7 need not reduce dilatational stresses 
so strongly, and so show no interesting microstructure. 

The key difference between the dynamics of ferroclas- 
tics and the dynamics of conventional systems is that 
the former is nonlocal. For example, the retraction of 
a streamer (or the vanishing of a domain) cannot be ef- 
fected by local changes; instead, a global change in the 
displacement field is,reauired. The nonlocality is seen ex- 
plicitly in theoriescJ'tZI where effective long-range inter- 
actions between the strains are generated on eliminating 
the displacement. And the relaxation stage is different 
also in that the changes from one snapshot to the next 
cannot be predicted by inspection. 

Figures 3(a), 3(b) and 3(c) show expanding regions of 
low-T phase, for stiffness parameter C = 10, 100 and 1000 
respectively. The 12-arm burst of part (a) is, we suggest, 
best viewed as a starburst disclination expanding into the 
high-T phase; it resembles strongly the "stellate twin" of 
Fig. 1 of Ref. U (though the latter is fully developed), 
consistent with our identification of cordierite as a soft 
system. The stiffer systems of parts (b) and (c) evolve 
more slowly; they also show branching at the outer edges, 
as discussed above. 

Figure 3(d) shows the system of 3(b) at a later time. 
There are several streamers, which retract very slowly 
with time; similar structures are seen in Fig. 14(b) of 
Ref. ^. The long timescale associated with this simple 
evolution (there is no change in the topology of the do- 
main walls) results from the nonlocal relaxation. 

Figures 3(e) and 3(f) show the system in 3(c) at a 
later time. The complicated structures in 3(e) are not 
observed in softer systems. In stiff systems they are tran- 
sitory, existing only immediately after completion of the 
phase transformation. The domain walls near the stri- 
ations and star embryos are relatively broad and they 
deviate from their optimal orientations; at later times, 
for the most part, they are replaced by well-defined walls 
as shown in the fully relaxed configuration of 3(f). The 
divided diamond-shaped structure seen at the top of Fig- 
ure 3(f) is seen in Fig. 9 (especially part (c)) of Ref. M, 
suggesting that Ta-N is a stiff system. 

VI. SUMMARY 

A Landau theory of proper hexagonal-orthorhombic 
ferroelastics, with the fewest possible number of pa- 
rameters, describes qualitatively many of the structures 
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seen in Ta4N. Mg-Cd alloys and Mg-cordicrite, and in 
trigonal-monoelinie lead orthovanadate. Our static and 
dynamical studies found structures not previously seen 
theoretically, including trapped high-temperature phase, 
12-armed bursts, streamers and striations. 

Since the theoretical understanding of structures in 
these systems seems now well advanced, we suggest that 
the next major theoretical task is to include noise and 
fluctuation effects in the time development. The effect of 
the inertial term also seems worth investigation. 



eh + + {e'h? = el, + {e',,f + {e'i.f , (23) 

with i ^ j. 

The isotropic case corresponds to ^2 = d'^; then I 
equals the ^2 invariant we used in the free energy. The 
^4 invariant is found from the difference, d'2 — —d'^. In 
keeping with our goal of a minimal theory, we consider 
the isotropic case and omit the d^ term in the free energy. 
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APPENDIX 

Equation (2) contains only one term quadratic in the 
strain gradients; the most general such expression in- 
cludes four additional invariants: 

+ ^3 [ei,i (e2,i + 63,2) + ei,2 (-62,2 - 63,1)] 



(i4 



(62,1-63,2) -(62,2 + 63,1)' 



+ ^5 (62,163,2 - 62,263,1) 



(20) 



The di term may be important more generally, but it 
can be neglected for the wall structures of interest to us 
here; it obviously plays no role for stiff systems where |ei| 
is small, and we find numerically that it is qualitatively 
unimportant also in soft systems. Given this, and that 
we want a minimal theory, wc omit the ^3 term also. The 

term integrates to a surface term, which vanishes for 
periodic boundary conditions. 

The remaining term (coefficient ^4) is closely related 
to the d2 term in Equation (2). To demonstrate this, we 
consider the sum 



1= ^4 [ell + (6;,i)' + (4in 

+ \d',[el2 + {e',,,f + {e'i,f] 



(21) 



of two non- negative invariants. The primes represent 
terms obtained by 27r/3 rotations about the hexagonal 
axis; for example, 63 ^ = \{—di + \/^d2)\{\/^e2 — 63). 
Each of the above six terms corresponds to the change 
in the order parameter in the direction normal to one 
of the six different domain- wall orientations discussed in 
Section III. The explicit expression for I is 



I=i^d'Meli + el2) + ^li + 4,2 



+2(62,163,2 + 62,263,1)] 

+ ^<[4l + e2,2 + 3(6L+62_2) 



-2(62,163^2 + 62,263,1)] • 



(22) 



Equivalent results are found on working with 62 rather 
than 63, for 
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COLOUR FIGURES 



FIG. 2. Samples of static configurations. The three differ- 
ent orthorhombic variants are represented as red, blue and 

green; the high-temperature phase is represented as black. 
Part (a) shows a very stiff system ((^ — 1000) at low temper- 
ature (r = —50); the microstructure is dominated by nested 
stars and (1,3,5,3) zero-disclination nodes. Part (b) shows an- 
other very stiff system ((^ — 1000), but at r = 0.8, just below 
the transition temperature {t — 1). The important features 
are the triangular and trapezoidal regions of trapped high-T 
phase and the (1,1,1,1,1,7) fan nodes. 



FIG. 3. Snapshots of dynamical simulations at different 
times for identical starting conditions and different values of 
the stiffness parameter, at r = —50 (well below the tran- 
sition temperature). Each frame has size 51.2 x 51.2; the 
colour scheme is as in Figure 2. Parts (a), (b) and (c), with 

= 10, 100 and 1000 respectively, show early-time structures 
(t=0.255) evolving from the perturbed high-T phase. The 
transformed region expands and then collides with images of 
itself because of the periodic boundary conditions; at later 
times the transformed phase occupies the whole frame. Part 
(d), a continuation of part (b) at time t = 8.191, shows sev- 
eral naxrow domains or "streamers". Parts (e) and (f) are 
continuations of part (c) at times t = 0.511 and t = 16.383 
respectively. Part (e) shows transitory structures associated 
with stiff systems, namely bi-domain striations and star em- 
bryos. Part (f) is a fully relaxed configuration. 
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